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ABSTRACT 

Aims. In this paper we investigate theoretical pulsation models for the d Scuti star 44 Tau. The star was monitored during several 
multisite campaigns which confirmed the presence of radial and nonradial oscillations. Moreover, its exceptionally low rotational 
velocity makes 44 Tau particulary interesting for an asteroseismic study. Due to the measured log g value of 3.6 ± 0. 1, main sequence 
and post-main sequence models have to be considered. 

Methods. We perform mode identification based on photometric and spectroscopic data. A nonadiabatic pulsation code is used to 
compute models that fit the identified modes. The influence of different opacity tables and element mixtures on the results is tested. 
Results. The observed frequencies of 44 Tau can be fitted in both the main sequence and the post-main sequence evolutionary stage. 
Post-main sequence models are preferable as they fulfill almost all observational constraints (fit of observed frequencies, position in 
the HRD and instability range). These models can be obtained with normal chemical composition which is in agreement with recent 
spectroscopic measurements. The efficiency of envelope convection (in the framework of the mixing-length theory) is predicted to 
be very low in 44 Tau. We show that the results are sensitive to the choice between the OPAL and OP opacities. While the pulsation 
models of 44 Tau computed with OP opacities are considerably too cool and too faint, the use of OPAL opacities results in models 
within the expected temperature and luminosity range. 
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1. Introduction 

6 Scuti stars are pulsating stars situated at the lower part of 
, the classical instability strip in the HR diagram. Their spec- 
tral types range from AO to F2 and they are in the hydrogen- 
core or hydrogen-shell burning stage. With masses between 1.5 
. and 2.5 Mq their periods of pulsation are between 30 min and 
■ 8 hours. In these stars pulsation is driven by the opacity mecha- 
I nism {k mechanism) acting in the Hell ionization zone. 

Among the 6 Scuti stars two subgroups are distinguished: 
' high-amplitude 6 Scuti stars (HADS) with peak-to-peak ampli- 
, tudes > 0.3 mag in V with mainly radial pulsation and low- 
amplitude 6 Scuti stars which oscillate predominantly non- 
radially. HADS stars show rather slow rotation in comparison 
with the average vsin/ of 96 km s"' of low amplitude 6 Scuti 
stars (Solano & Fernlev ,1997) . 44 Tau (HD 26322, spectral type: 
F2 IV) exhibits characteristics of both groups: radial as well as 
nonradial modes are excited to observable amplitudes and its 
projected rotational velocity is very small (vsin/ = 2±lkms"'). 
As discussed by Antoci et al. ( 120071 ) (hereafter ABR07) this in- 
dicates either pole-on view and/or intrinsic slow rotation. Zima 
et al. (120071) recently confirmed that 44 Tau is an intrinsic slow 
rotator This makes this star particulary interesting for astero- 
seismology. 

Previous theoretical studies on this star have been carried 
out by Civelek et al. ('2001') and Kirbiyik et al. ( 12003) . By means 
of fitting adiabatic oscillation frequencies to the observed data 
they conclude that 44 Tau is a 1.89 M© star at the hydrogen- 
shell burning stage. Recently, Garrido et al. (i2007j concluded 
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that the pulsations of 44 Tau could be explained by a main se- 
quence model with a mass of 1.94 M©. 

This paper extends the work of Lenz et al. ilQQl ) and 
presents a comprehensive seismic study based on the detailed 
frequency analysis of ABR07 and on additional photometric data 
obtained in 2004/05. 

The paper is structured into the following parts: Section |2] 
summarizes the results of available photometric and spectro- 
scopic data. Section |3] is devoted to the theoretical modelling of 
the frequencies of 44 Tau. After the determination of the spher- 
ical degrees of the observed frequencies we compute models to 
fit the detected frequencies and the observed frequency range of 
excited modes. Finally, the effect of using different opacity ta- 
bles and element mixtures is studied. In Section |4] we present 
our conclusions. 

2. Observations 

2.1. Photometry 

44 Tau was the main target of several Delta Scuti Network 
(DSN) campaigns. Photometric data obtained during three cam- 
paigns between 2000 and 2003 were analysed in detail by 
ABR07. 29 significant frequencies of which 13 are independent 
were detected. 

In the winter season 2004/05 an additional observing run 
was dedicated to this star. The new photometric data confirmed 
the previous results. Although the noise for the total data set is 
much lower only a few additional frequencies could be detected. 
However, these frequencies are combinations of already detected 
frequencies and no additional independent frequency was found. 
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Table 1. Averaged observed amplitude ratios and phase differences from 2000/1 to 2004/5. Only frequencies detected in both 
photometric and spectroscopic data are listed. In the last two columns the predicted spherical degree from our mode identification 
(see Section [34] | and the azimuthal number m from Zima et al. (2007) are given. For fg no reliable results could be derived due to 
the influence of the close and more dominant frequency f(,. 
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Some modes in 44 Tau show strong amplitude variability. 
This is not unusual and is a rather common phenomenon in 
6 Scuti stars (Poretti 120001 Breger 2000 ). There are several possi- 
ble explanations such as, for example, beating of close frequen- 
cies or resonance effects. In Fig.[T]the frequency spectrum of the 
2004/05 data is compared to the results from the 2000/01 data 
set. The presence of amplitude variations is clearly visible. 
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Fig. 1. Comparison of Stromgren v amplitudes from the observ- 
ing season 2000/01 with the data from 2004/05. The error bars 
are smaller than the size of the given points. The frequencies f2, 
f3, f4 and f6 show strong amplitude variation. 



2.2. Spectroscopy 

In 2004 a multisite spectroscopic campaign was carried out. 
Zima et al. ( I2007I I detected 16 frequencies of which 4 are com- 
binations. 12 independent frequencies were found in the spec- 
troscopic as well as in the photometric data. Using the moment 
method they constrained the azimuthal order of 8 modes, fi, f4, 
fs, f7, fg and fll were identified as axisymmetric modes whereas 
f2 and f3 were found to be prograde modes with m - 1. The 
corresponding frequency values are listed in Table [1] 

Only two modes were determined to be non-axisymmetric, 
therefore only a rough determination of the inclination angle / 



was possible. Zima et al. (120071 ) estimated / to range between 
35° and 85°. Since v sin / is derived to be 2 + 1 km s"' this would 
correspond to an equatorial rotational velocity between 1 and 
5 km s"'. Note that the newly derived value for the projected 
rotational velocity is slightly lower than the results of previous 
studies by Solano & Fernley d 19971 1 who measured v sin / to be 
6.8 ± 1.2 km s"'. 

The analysis of the photospheric abundances of 44 Tau by 
Zima et al. (2007) shows a similar composition as for the sun. 
The effective temperature is estimated to be 7000 ± 200 K and 
log^ is determined to be 3.6 ± 0.1 consistently from line profile 
and equivalent width measurements. 

3. Seismic l\/lodelling 

3. 1 . Fundamental Parameters of 44 Tau 

Based on the measured HIPPARCOS parallax of 
16.72 ± 0.93 mas, a luminosity of logL/Lo = 1.340 ± 0.065 
was derived. The effective temperature was estimated from 
Stromgren and Geneva photometry and from the Vienna New 
Model Grid of Stellar Atmospheres (NEMO, Nendwich 
et al. 120041 Heiter et al. .20U2l l. The derived value is 
logTefF = 3.839 + 0.007 (6900 K + 100 K). This result is 
in agreement with the value inferred from spectroscopy. 

The spectroscopically measured log g value of 3.6 + 0.1 does 
not allow an unambiguous determination of the evolutionary sta- 
tus of 44 Tau. This can be clearly seen in a \ogg vs. log T^ff plot 
such as Fig. 5 in Pamyatnykh (2000). Consequently, both main 
sequence and post-main sequence models have to be considered. 

Since the intrinsic rotational velocity of 44 Tau is very low, 
we consider only nonrotating models in our study. 

3.2. Tools 

The models in this paper were computed with an improved ver- 
sion of the codes described by Breger & Pamyatnykh (19981. 
The OPAL opacity tables (Iglesias & Rogers, |19961 l supple- 
mented with the low-temperature data of Alexander & Ferguson 
( 119941 ) were used. If not stated otherwise, we made use of the el- 
ement mixture from Grevesse & Noels (1993 hereafter GN93) 
and the last version of the OPAL equation of state (see Rogers 
et al. 1 19961 version OPAL EOS 2005). Nuclear reaction rates are 
taken from Bahcall and Pinsonneault ( I1995I I. A linear nonadia- 
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Fig. 2. The Hertzsprung-Russell diagram showing the location of 
models that fit radial fundamental and first overtone modes to the 
observed frequencies for the main sequence (upper panel) and 
post-main sequence case (lower panel). The asterisk indicates 
the center of the photometric error box. Grey points represent 
models with Z values between 0.019-0.04 (MS) or 0.015-0.025 
(post-MS) and an overshooting parameter q-qv ranging from to 
0.4. Representative evolutionary tracks for the models listed in 
Table |2] are plotted. 



Fig. 3. Petersen diagrams for a main sequence (upper panel) and 
a post-main sequence model (lower panel). The asterisk corre- 
sponds to the observed values of 44 Tau. The uncertainties in 
period (~ 1.4 x 10"^ cd"') and period ratio (~ 2 x 10"^) are 
smaller than the symbol size. 



scale heights. As a representative example, the Petersen diagram 
for two selected models is shown in Fig. |3] The corresponding 
model parameters are listed in Table |2] 



batic analysis of low-degree oscillations was performed using a 
modified version of Dziembowski's code (1977). 

3.3. Radial Modes 

In ABR07 the radial character of the frequencies fi and fs was 
confirmed. The fact that radial fundamental and first overtone are 
observed can be used to reduce the number of possible models 
significantly. To find appropriate models for 44 Tau it is useful 
to utilize Petersen diagrams (Petersen & Jorgensen |1972| l. In the 
Petersen diagram the period ratio of the radial first overtone to 
the fundamental mode is plotted as a function of the period of the 
fundamental mode. These diagrams can be used to determine the 
correct mass of the model which fits the radial fundamental fre- 
quency and the ratio F/IH (fi/f5=0.7698) to the observed values. 
The period ratio is sensitive to metallicity and rotation as shown 
by Suarez et al. (2006). In Section [378] we will demonstrate the 
effect of different opacity tables and the new solar element mix- 
ture on the period ratio. 

We have computed stellar models for a large set of different 
input parameters such as chemical composition, mixing-length 
parameter, omlt, and the parameter for overshooting from the 
convective core, aov Our results show that main sequence and 
post-main sequence models which fit the radial fundamental and 
first overtone modes are clearly separated in the HR diagram 
(see Fig. |2]i. For each case an evolutionary track of a selected 
model is delineated. Post-main sequence models are located in- 
side the photometric error box, whereas main sequence models 
are predicted at significantly lower effective temperatures in the 
HR diagram. The given family of models was computed using 
a wide range of parameters, such as Z between 0.015 and 0.04 
and with an overshooting parameter of up to 0.4 pressure 



3.4. Identification of Nonradial Modes 

In the previous section we have noted that the observation of two 
radial modes puts strong constraints on the models. For complete 
seismic modelling it is necessary to obtain an identification for a 
sufficient number of observed nonradial pulsation modes. 

Mode identification is often complicated by the effects of ro- 
tation. The exceptionally small rotational velocity of 44 Tau al- 
lows to reduce these problems because the rotational splitting is 
very small. There are several techniques that either make use of 
photometric data in two filters or combine them with spectro- 
scopic data to infer the spherical degrees of the modes. The lat- 
ter approach has been successfully applied in the case of FG Vir 
(Daszynska-Daszkiewicz et al. 12005). The azimuthal order of 
pulsation modes can only be determined from spectroscopy. 

To obtain reliable photometric mode identification we com- 
puted average amplitude ratios and phase differences in the 
Stromgren filters v and y using weights for each observing sea- 
son to account for the different amount and quality of the data 
(see formulae in the Appendix of Breger et al. I1999I I. The ob- 
servational uncertainties were derived by error propagation from 
the uncertainties of the annual solutions. These results are given 
in Table □ 

For 6 Scuti stars theoretical mode positions in the ampli- 
tude ratio vs. phase difference diagram are sensitive to the ef- 
ficiency of convection in the stellar envelope. Montalban & 
Dupret (12007) showed that different treatments of convection af- 
fect the phase differences and amplitude ratios as well. In our 
computations, we use the mixing-length theory of convection 
and assume the frozen convective flux approximation for the 
pulsation-convection interaction. This assumption may not be 
adequate for cold 6 Scuti stars and result in large uncertainties in 
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Table 2. Physical parameters of two representative models used in this paper. If no indication is given in figures the corresponding 
plots refer to the models in this table. 



X 



ttMLT ffov M/Mo log Teff log L log g Age [Myr] 



MS 0.70 0.03 0.0 0.25 2.010 3.8077 1.2385 3.6848 1150 

post-MS 0.70 0.02 0.0 0.0 1.875 3.8422 1.3601 3.6712 1120 
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Fig. 5. Observed mode positions with error bars. For better clarity modes with small and large error bars are shown in separate 
diagrams. In both diagrams the predicted mode positions (for modes within the observed frequency range) refer to the case of 
inefficient convection, ckmlt - 0.0. The diagrams plotted for a stellar model with 1.85 M©, X=0.70, Z=0.02 and q-qv = 0.0. 
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Fig. 4. Diagnostic diagram used for mode identification. 
Comparison of theoretical results for qtmlt - 0.5 and umlt = 0.0 
to the observational positions. The asterisks correspond to the 
amplitude ratios and phase differences of the observed frequen- 
cies. Theoretical mode positions are given as points connected 
with lines. Only predicted modes within the observed frequency 
range are plotted. 



the mode identification. However, Montalban & Dupret ( I2007I I 
find that for frequencies close to that of the radial fundamental 
mode the convection-pulsation treatment does not influence the 
determination of mode degree, {. 

We tested diff'erent values of the mixing-length parameter 
^MLT for a stellar model with standard chemical composition 
and a mass of 1.85 Mq. The results of these computations are 
shown in Fig.|4]and|5]. The theoretical mode positions were cal- 
culated according to Daszynska-Daszkiewicz et al. (I2003I I using 
photometric data in two passbands. 



As can be seen in Fig. |4] for q-mlt - 0.5 (or higher) there 
is no reliable agreement between observed and theoretical mode 
positions. Instead, we obtain a good consistency for inefficient 
convection (q-mlt ^ 0.2). In the two diagrams in Fig.|5]the posi- 
tions of observed modes are compared to the theoretical predic- 
tions for a model with inefficient convection. We can definitively 
identify the frequencies fi and fg as {=0 modes, fi, h, f4 and fe 
as £-1 modes, fj can be identified as an ^ = 2 mode. The ob- 
servational formal errors for the remaining frequencies are too 
large to draw a conclusion. However, an identification as ^ = 2 
is probable for some of ffie modes. Still, it is necessary to con- 
firm their mode degree by anoffier method. 

Since radial velocity data from 2004 are available (Zima et 
al. I2007I I. it is possible to combine them wiffi photometric data 
from the Stromgren v and y passbands using ffie method by 
Daszynska-Daszkiewicz et al. ( |2003l|2005] l. This method allows 
to perform a^^ test which indicates the most probable spherical 
degree for a mode. The results agree with those from the purely 
photometric mode identification. Additionally this method iden- 
tifies fg and fio as ^ = 2 modes with a high confidence level. 
Table[T]summarizes the results of our mode identification. 

As mentioned before, we rely on the assumption of frozen 
convective flux, which may not be a realistic assumption for 
cold S Scuti stars. A comparison with mode identification re- 
sults using a time-dependent convection approach is thus de- 
sirable. Recently Garrido et al. (I2007I I computed diagnostic di- 
agrams to identify the pulsation modes of 44 Tau using the 
time-dependent convection approach described by Dupret et al. 
(I2004I I and Grigahcene et al. ( I2005I I. Even though using diff'er- 
ent treatments of convection ffieir results are in agreement with 
our results. Within the same context but for another 6 Scuti star 
(FG Vir) Daszynska-Daszkiewicz et al. f2005i l have compared 
the results for the complex parameter / and mode instability 
from the frozen convective flux assumption and ffie nonlocal 
time-dependent mixing-length theory (Gough |1977a,b| l. As for 
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44 Tau, envelope convection is inefficient in FG Vir and the re- 
sults from both treatments are similar. Consequently, our simple 
approach with a small mixing-length parameter (ffMLX 5 0.2) 
may also provide reliable results in the case of 44 Tau. 

By comparing our mode identification results with the fre- 
quency spectrum in Fig. [T] an interesting conclusion can be 
drawn: the strongest amplitude variability is shown by ^ = 1 
modes, whereas the two radial modes are constant in their am- 
plitudes. This fact may be important to explain the cause for the 
observed amplitude variability in 44 Tau. 

Two of the {-1 modes, f3 and fe, are separated by only 
0.44 cd"' whereas the frequency distance between the other 
{=1 modes fa and or fg and f4 is around 2 cd~'. This pat- 
tern indicates that we observe an avoided crossing phenomenon 
(Aizenman et al.| 19771), which is a coupling between two modes 
of different type, viz. an acoustic (p) and a gravity (g) mode. If an 
avoided crossing is indeed observed it represents an additional 
strong constraint to the model (see Dziembowski & Pamyatnykh 
11991 ). In Fig. |6]the evolutionary change of the frequencies is 
shown. The given box restricts the uncertainties of the measured 
effective temperature and the observed frequency range. Within 
this box an avoided crossing between gi and g2 takes place (the 
gi mode behaves here like an acoustic p-mode due to the previ- 
ous avoided crossing with the pi mode). For the main sequence 
model which fits the radial modes (vertical line) we predict two 
avoided crossing pairs for { = I within the observed frequency 
range. One of them fits the observed frequencies as we will see 
later. 



3.5. Instability Ranges 

Successful seismic modelling also involves the correct predic- 
tion of instability of the detected modes. Instability reflects the 
properties of the envelope, where mode driving occurs. 

Instability ranges for two representative models (on the main 
sequence and in the post-main sequence stage) which fit the ra- 
dial as well as i-l modes are shown in Fig. |7] In this diagram 
the normalized growth-rate as defined by Stellingwerf (I1978I I, 

rR rR 

^ — Jq (dW/dr) dr / |dW/dr| dr (where W is the usual 
work integral), is plotted against frequency, r] varies between - 
1 (full damping) and -i-l (full driving). It can be seen that for 
the main sequence model the instability range (rj > 0) covers 
the observed frequency range completely but extends to consid- 
erably lower frequencies than observed. For the post-main se- 
quence model the instability range extends to higher frequencies 
than actually observed. The theoretical lower boundary for un- 
stable frequencies fits the observational limit much better than 
in the main sequence case but the lowest observed frequency fio 
is predicted to be stable. With exception of a significant effect 
of convective efficiency on the high frequency boundary of the 
instability range, the results are not very sensitive to changes of 
the input parameters. 

Generally, the value of the highest unstable frequency is very 
sensitive to the efficiency of convection and to the convection- 
pulsation treatment. Since our code relies on the simplest frozen 
convection assumption, it is not suitable to predict the upper 
limit of the frequency instability range. If future observations re- 
veal additional independent frequencies that extend the observed 
frequency range, further constraints on the evolutionary status 
and on the efficiency of convection in 44 Tau may be provided. 
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Fig. 7. Comparison of the observed and theoretical frequency 
ranges for main sequence (upper panel) and post-main sequence 
(lower panel) models. The normalized growth-rate, rj, is positive 
for unstable modes. The vertical lines represent the observed fre- 
quencies (with longer lines for { - modes). 



3.6. Main Sequence Models 

As discussed in Section [J!2l it is not possible to fit the observed 
period ratio Pi/Po and the measured fundamental parameters si- 
multaneously for models in the hydrogen-core burning stage. 
The predicted main sequence models for 44 Tau are generally 
too cool (see Fig.|2]i. However, the photometric calibrations may 
have larger uncertainties than adopted in Fig. |2l Therefore, we 
also consider MS models for 44 Tau. We computed a grid of 
models by varying the chemical composition and the overshoot- 
ing efficiency. Not for all models the radial fundamental and first 
overtone mode could be fitted on the main sequence. Using stan- 
dard chemical composition (X=0.70, Z=0.02) and assuming in- 
efficient convection, q-mlt = 0.0, it is not possible to obtain any 
main sequence model for q-qv < 0.25. For aov = 0.25 the model 
is located exacdy at the TAMS. For Z=0.019 and 0.021 main 
sequence models are possible for ofov > 0.28 and 0.23, respec- 
tively. Thus, the presence of effective convective overshooting is 
essential to model 44 Tau in the stage of hydrogen burning in the 
core. 

The evolutionary track for the main sequence model with the 
best fit of the identified observed frequencies is delineated in the 
upper panel of Fig.|2] This 2.01 Mq model was computed using 
X=0.70, Z=0.03, Q'MLT=0.0, Q'ov=0.25. Additional information 
on the physical parameters of the model is given in Table|2] The 
predicted frequency spectrum (see Fig. [8] upper panel) shows 
good agreement between predicted and observed frequencies of 
{ - Q and { = I modes. A good fit of the { = 2 modes could not 
be obtained. 

3. 7. Post-Main Sequence Models 

In contrast to main sequence models it is not difficult to obtain 
models within the photometric error box in the post-main se- 
quence case (Fig.|2] lower panel). The post-main sequence evo- 
lution across the observed error box is only 5-9 times faster than 
the main sequence evolution. Therefore, the probability to ob- 
serve a star in the H-shell burning stage is not small. For the rep- 
resentative model given in Fig.|2l the input parameters X=0.70, 
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Fig. 6. Behaviour of the frequencies of ^ = 1 modes (left panel) and { - 2 modes (right panel) in the course of the main sequence 
evolution from the ZAMS (leftmost points) to the TAMS (rightmost points). Larger points represent unstable modes. The vertical 
line represents the effective temperature of the model with the best fit to the observed frequencies. The box indicates the ranges for 
observed frequencies and for effective temperature determined from photometry. The parameters of the model are given in Table |2l 
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Fig. 9. Kinetic energy density of pulsation along the stellar ra- 
dius for three £ = 1 modes. The modes at 9.21 and 11.65 cd"' 
are trapped, whereas the mode at 10.25 cd"' is not. The energy 
contribution from the outer envelope is significant for the two 
trapped modes. 



Z=0.02, aMLT^O.O, aov=0.0 and an initial mass of 1.875 
were chosen. The corresponding frequency spectrum for this 
model which fits the observed radial modes can be seen in Fig. [8] 
(lower panel). 

The frequency spectrum of nonradial modes is much denser 
than that of main sequence models. In fact, for post-main se- 
quence models many more frequencies are predicted than ob- 
served. An extreme example is 4 CVn, one of the most evolved 
6 Scuti stars, for which approximately 500 modes are predicted 
(Breger & Pamvatnvkh |20021 l. There certainly exist some rules 
which select only a few modes to be excited to observable am- 
plitudes. Mode trapping in the acoustic cavity in the stellar en- 
velope can be one of these rules, as discussed by Dziembowski 
& Krolikowska ( fT990] l. 

Trapped modes can be considered as nonradial counterparts 
of radial modes. Most of the energy of trapped modes is confined 



in the envelope, where the density is low. Therefore, the total ki- 
netic energy of these modes is much smaller than that of modes 
which are not trapped. Consequently, it may be easier to excite 
trapped modes to observable amplitudes. In Fig. |9] the distribu- 
tion of the kinetic energy density of pulsation along the stellar 
radius for two trapped {=1 modes and for a mode which is not 
trapped is shown. 

Mode trapping can also explain the tendency of detected fre- 
quencies to build groups around radial frequencies as can be seen 
in Fig. [To] This diagram shows the increment of the exponential 
amplitude growth, y - -Im(w), where a> is the complex eigen- 
frequency. y is proportional to the work integral and inversely 
proportional to the kinetic energy of a given mode. Just the ki- 
netic energy is responsible for its nonmonotonic behavior. The 
increment y is significantly higher for trapped modes. 

We can also see, that trapping of ^ = 2 modes is rather weak 
compared to the £ = I and { - 3 case (see Fig. [TOl l. The rea- 



Lenz et al.: An asteroseismic study of the S Scuti star 44 Tau 



7 



l.b r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 

1.0 - I III I I 1 1 I II I 

0.5 - 




-4.0 - * trapped ( = 1 



-4.5 - -o-( = 2 

-5.0 - h"-"^ I - 

-5 5 [ I I I I I I I I I 1 I I I I I 1 I 1 I I I I 

6 7 8 9 10 11 12 13 14 15 16 17 

Frequency [cd"'] 

Fig. 10. Increment of the exponential amplitude growth, y, for 
unstable modes of ^ < 3 for the post-MS model shown in Fig. 6. 
The trapped {=1 modes are marked with black triangles. 



son for this becomes clear when the run of the Brunt- Vaisala 
frequency, A^, and the Lamb frequency, L, is examined in detail. 
These characteristic frequencies determine the regions with the 
oscillatory behaviour in the star as shown in Fig.[TT] Similar di- 
agrams can be found, for example, in Unno et al. 1979 (see Fig. 
15.1 there). The propagation zone for p-modes is defined by the 
conditions > and > where v (= Re(w)/(2;7r) ) denotes 
the mode frequency, g-modes are confined mainly in stellar re- 
gions where < A^^ and < L^. Between these two cavities 
there is a so-called evanescent zone with exponentially decreas- 
ing amplitude behaviour While the Brunt-Vaisala frequency is 
independent of the spherical degree, the Lamb frequency incor- 
porates the factor + 1). As can be seen in Fig.[TTl the mode 
cavities of p- and g-modes are poorly separated for £ - 2 modes 
in the observed frequency range (note the very small zone be- 
tween N and L2, in contrast to the more extended zones between 
Li and N for { = I and between N and L3 for { = 3). Therefore, 
no effective mode trapping occurs for f - 2 modes. 

This might explain why { = 2 modes are observed with rela- 
tively small amplitudes only in 44 Tau as well as in other 6 Scuti 
stars (such as FG Vir and 4 CVn). 

3.8. Effect of Opacities and Element Mixture 

As rotational effects can be neglected due to the low measured 
rotational velocity of about 3 + 2 km s"', 44 Tau is a good target 
to test the effects of different input data on the Petersen diagram. 

The results in the preceding sections were obtained with the 
OPAL opacity tables. Another, completely independent set of 
opacity tables is available from the Opacity Project (OP, see 
Seaton 2006 ). The use of different opacity tables slightly influ- 
ences the location and size of the opacity bumps inside stellar 
models and therefore influences the driving of pulsations. Due 
to differences in opacities we expect also differences in the stel- 
lar structure and in the predicted frequencies of pulsations. 

We compared the preceding results with those obtained by 
using the OP tables. Using the same model parameters and ele- 
ment mixture (Grevesse & Noels et al. 1993 ) as for the standard 
OPAL post-main sequence model, we computed evolutionary 
tracks and corresponding frequencies. In the Petersen diagram 
(Fig.[T2li the difference of the period ratio between models com- 
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Fig. 11. Run of the Brunt-Vaisala and Lamb-frequencies for 
{ < 3 inside the post-MS model of 44 Tau. The two horizon- 
tal lines mark the observed frequency range (v between 6.3 and 
12.7 cd"'). The extent of the evanescent zone is small for £ = 2 
resulting in only weak trapping of f = 2 modes. 

puted with the two opacity sets is shown. The period ratios for 
the OP case are higher than those obtained with OPAL tables. 
Hence for models computed with OP opacities a lower mass is 
predicted for 44 Tau. However, as can be seen in Fig. [T3j seis- 
mic models computed with the OP opacities are clearly too cool 
and too faint. We do not consider main sequence models here, 
because their effective temperatures are already very low using 
OPAL tables (Fig.|2l upper panel). 



0.778 



0.776 - 




-1.4 -1.3 -1.2 -1.1 -1.0 -0.9 -0.8 -0.7 



Fig. 12. Comparison of the results for OPAL and OP opacities 
in the Petersen diagram in the case of a nonrotating post-main 
sequence model with standard chemical composition, inefficient 
convection (ffMLT = 0.2) and assuming no overshooting from the 
convective core. The mass for the OP model was reduced to fit 
the observed period ratio of 44 Tau. 

Recently, Asplund et al. (I2004I I found evidence for a revi- 
sion of the solar chemical composition (see also Asplund et al. 
2005). Their results indicate that the solar abundances of carbon, 
nitrogen, oxygen and other heavy elements are lower than pre- 
viously estimated. We also tested the influence of the new solar 
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Fig. 13. Comparison of the evolutionary tracks of the models 
given in Fig. [12] 




Fig. 15. HR diagram showing the evolutionary tracks of the mod- 
els presented in Fig. [14] 



element mixture (hereby referred to as A04) on the OPAL mod- 
els of 44 Tau. In the Petersen diagram in Fig. [14] we compare our 
standard post-MS model with models computed with A04 mix- 
ture and using the new solar values X=0.74, Z=0.012. As can 
be seen, only a small change in mass has to be applied to fit the 
period ratio. However, during the main sequence evolution, the 
element mixture itself has a significant effect on the period ratio 
as well. In the HR diagram (Fig. [TSj both models can be fitted 
within or close to the photometric error box. Furthermore, we 
obtain the same mode identification as for the models computed 
with the GN93 mixture. The instability range is only slightly 
shifted to lower frequencies, still predicting frequency fio not to 
be excited. It is possible to fit the observed frequencies at al- 
most the same accuracy as for the GN93 mixture. Additional 
tests with main sequence models provided similar conclusions 
as in the case of post-MS models. 



0.778 



0.776 - 




log P. 



Fig. 14. Comparison of period ratios for GN93 and A04 element 
mixtures for post-main sequence models obtained with OPAL 
opacities. For the models computed with A04 the new solar 
values X-QJA, Z=0.012 were used. Just as in Fig. [12] a non- 
rotating post-main sequence model with inefficient convection 
(»MLT = 0.2) and no overshooting from the convective core was 
assumed. 



It is important to note that the effects of both opacities and el- 
ement mixture on the period ratio are comparable or even larger 
than rotational effects. 

4. Conclusions 

We presented a comprehensive asteroseismic study for the 
slowly rotating 6 Scuti star 44 Tau. The spectroscopically deter- 
mined log g value of 3 .6 + 0. 1 allows no definitive determination 
of the evolutionary stage, so the main sequence and post-main 
sequence case were considered. 

Our mode degree identification based on photometric and 
spectroscopic data determines fi and fs as radial fundamental 
and first overtone modes. This confirms the results of previous 
studies. Four frequencies, f2, f3, f4 and fe, were identied as ^ = 1 
modes, and three frequencies, f?, fg and fio, as f = 2 modes. We 
used the simplest model of convection (mixing-length theory) 
and the simplest convection-pulsation treatment (frozen convec- 
tive flux assumption) which may not be adequate for cold 6 Scuti 
stars, as discussed by Montalban & Dupret (120071 1 . However, our 
results are in good agreement with those obtained by Garrido et 
al ( 120071 ) using time-dependent convection. 

The identification of the radial fundamental and first over- 
tone frequency poses an important constraint to the models as it 
reduces the number of possible models significantly. Moreover, 
it is likely that an avoided crossing of nonradial modes in 44 Tau 
is observed. This fact puts additional strong constraints on the 
parameter space for appropriate models. For both considered 
evolutionary stages it is possible to fit the observed radial period 
ratio. 

If standard chemical composition and no overshooting from 
the convective core are applied, only post-main sequence mod- 
els of 44 Tau are located within the photometric error box in the 
HR diagram. Main sequence models can only be obtained by 
assuming convective overshooting. Examining such hydrogen- 
core burning models within a wide parameter range we find that 
all these models are significantly too cool. This fact is a major 
problem for main-sequence models. Nevertheless, it is remark- 
able that all four identified { - I modes can be fitted in both 
evolutionary stages. 

The predicted range of unstable modes is in good agreement 
with the observed frequencies in both cases. For main sequence 
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models the instability range extends to lower frequencies than 
observed while for post-main sequence models higher frequen- 
cies are found to be unstable. The lowest frequency fio is pre- 
dicted to be stable in the post-main sequence models, but it is 
close to the theoretical instability region. Our results on the in- 
stabiUty ranges can be used to put a further constraint on the 
determination of the evolutionary stage if new frequencies are 
detected that extend the frequency range. Additional frequen- 
cies beyond the highest observed frequency may also help to 
constrain the efficiency of convection more precisely. 

Considering the good fit of the observed frequencies and 
physical parameters such as eff^ective temperature, we find that 
the post-main sequence stage is preferable to explain the pulsa- 
tion spectrum of 44 Tau. 

In this investigation, we also tested the effect of opaci- 
ties and element mixtures on the radial modes. The choice be- 
tween OPAL and OP tables causes a significant difiference in the 
Petersen diagram. In the case of 44 Tau, the OPAL opacities are 
preferable for the asteroseismic analysis, because with the OP 
data it is not possible to fit the radial modes of 44 Tau within 
or close to the observational error box in the HR diagram. The 
use of the new solar element mixture results in a slightly cooler 
post-main sequence model. We do not find indications in favor 
of the old or the new element mixture. However, it is important 
to note that the influence of different opacity data on the period 
ratio should also be considered when using the Petersen diagram 
as a diagnostic tool. 
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